Shot noise in resonant tunneling through an interacting quantum dot with intradot 

spin-flip scattering* 
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In this paper, we present theoretical investigation of the zero- frequency shot noise spectra in 
electron tunneling through an interacting quantum dot connected to two ferromagnetic leads with 
possibility of spin-flip scattering between the two spin states by means of the recently developed 
bias- voltage and temperature dependent quantum rate equations. For this purpose, a generalization 
of the traditional generation-recombination approach is made for properly taking into account the 
coherent superposition of electronic states, i.e., the nondiagonal density matrix elements. Our 
numerical calculations flnd that the Fano factor increases with increasing the polarization of the 
two leads, but decreases with increasing the intradot spin-flip scattering. 



I. INTRODUCTION 



Recent rapid development of spintronic and single 
electron devices has resulted in intensive studies of 
spin-related phenomena in various resonant tunneling 
structuresi^ One of the devices of interest is a quan- 
tum dot (QD), sandwiched between two ferromagnetic 
electrodes i^i^ This kind of spin-related single-electron 
devices suffers inevitably intrinsic relaxations (deco- 
herence) due to the spin-orbital interaction^ or the 
hyperfine- mediated spin-flip transition^ The effect of the 
intrinsic decoherence on the quantum transport proper- 
ties in these devices still deserves more investigation. As 
is well known, measurement of shot noise, which is a 
consequence of the discrete nature of the electric charge, 
can provide information about the microscopic interac- 
tions and the statistics of particles, which is not available 
through the conductance measurement alonei^ There- 
fore, it is remarkably desirable to expand our knowledge 
about the shot noise spectra in such a spin-related single- 
electron tunneling device. 

For metallic QD without intrinsic decoherence, the 
zero-frequency Fano factor, which quantifies correlations 
with respect to the uncorrelated Poissonian noise, was 
analyzed with the classical rate equation approach (clas- 
sical shot noise) in Refs. pl]irill ll2] . However, in order to 
account for the quantum effects during sequential tunnel- 
ing processes, such as coherence and superposition of the 
wave functions and various intrinsic interactions, a more 
elaborate quantum description is required. The quantum 
rate equation approach is a suitable tool for accomplish- 
ing this task&iiiii^ In particular, some of the authors, 
Dong, Cui, and Lei, recently established the quantum 
rate equations for coherent tunneling through the cou- 
pled quantum systems with a bias-voltage and tempera- 
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ture dependent version^ Consequently, these equations 
allow us to study the quantum tunneling in moderately 
small bias-voltage region, where the quantum coherence 
plays more prominent role. Therefore, the purpose of the 
present paper is to develop an approach for calculation 
of the coherent effect, due to the intrinsic spin scattering, 
on the zero-frequency shot noise spectra in an interact- 
ing QD weakly connected with two ferromagnetic leads 
based on the quantum rate equations. 

The rest of this paper is organized as follows. In section 
II, we introduce the Hamiltonian for the interacting QD 
with intradot spin scattering, and then give the quantum 
rate equations. In section III, we describe our methodol- 
ogy for evaluating the quantum shot noise. In section IV, 
we present the numerical results and discussions. Finally, 
a brief summary is presented in section V. 



II. SINGLE QUANTUM DOT: MODEL AND 
QUANTUM RATE EQUATIONS 

The system that we study is a single QD with an ar- 
bitrary intradot Coulomb interaction U connected with 
two ferromagnetic leads. In this paper, we assume that 
the tunneling coupling between the dot and the leads is 
weak enough to guarantee no Kondo effect and that the 
QD is in the Coulomb blockade regime. For simplicity, we 
model the intrinsic spin relaxation with a phenomenologi- 
cal spin-flip term R and assume that this spin-flip process 
just happen in the dot (the spin conserving tunneling). 

Moreover, we assume that the temperature is small 
enough to see the effects due to the discrete charging and 
discrete structure of the energy levels, i.e. ksT <C U,A, 
where, A is the energy spacing between orbital levels. 
Each of the two leads is separately in thermal equilib- 
rium with the chemical potential /i,j, which is assumed 
to be zero in the equilibrium condition and is taken as the 
energy reference throughout the paper. In the noncqui- 
librium case, the chemical potentials of the leads differ 
by the applied bias V. We are interested in the regime 



eV <C A, where only one dot level (ed) contributes to 
the transport. Here we neglect Zeeman splitting of the 
level due to weak magnetic fields B {gusB < ksT), 
which means that both the spin-up and spin-down trans- 
ports through the dot go through the same orbital level 
Cd- Therefore, the Hamiltonian of resonant tunneling 
through a single QD can be written as>i^ 



H ^ Y^ e^fca4fe<,c^,fe^ +^dY^ cl„c^„ + R{c\ 
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where c j.^ {cr^ka) and cjj^ [cdc,) are the creation (annihi- 
lation) operators for electrons with momentum fc, spin-a 
and energy €r,ka in the lead 77 (= L, R) and for a spin-a 

electron on the QD, respectively. Uda = Cda'^da '^^ ^^'^ ^'"~ 
cupation operator in the QD. The fourth term describes 
the Coulomb interaction among electrons on the QD. The 
fifth term represents the tunneling coupling between the 
QD and the reservoirs. We assume that the coupling 
strength Vr^a is spin-dependent being capable of describ- 
ing the ferromagnetic leads. 

Under the assumption of weak coupling between the 
QD and the leads, and applying the wide band limit in 
the two leads, electronic transport through this system 
in sequential regime can be described by the bias- voltage 
and temperature dependent quantum rate equations for 
the dynamical evolution of the density matrix elementsii^ 



Poo — 2^(r^crPo-o- - ^rjaPOo), 
r}(7 
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pdd = Y,KiP^^ ^Y^K^p'^i^Y.^^vi +^vOP'i<i' (2d) 
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(cr =1', J, stands for electron spin and a is the spin oppo- 
site to a). The statistical expectations of the diagonal 
elements of the density matrix, pa (i = {0, cr, d}), give 
the occupation probabilities of the resonant level in the 
QD as follows: poo denotes the occupation probability 
that central region is empty, p^a means that the QD 
is singly occupied by a spin-a electron, and pdd stands 
for the double occupation by two electrons with differ- 
ent spins. Note that they must satisfy the normaliza- 
tion relation poo + Pdd + J^a Pec — 1- The non-diagonal 



elements p^a describe the coherent superposition of dif- 
ferent spin states. These temperature-dependent tunnel- 
ing rates are defined as F^^ = ^jja-f^i^d) and r^^. = 
Vriaf^i^d + U), where r,,cr are the tunneling constants, 
/+(cj) = {1 -^ ei^-t^v)/Tyi is |;i^g Fermi distribution 
function of the rj lead and /^(w) = 1 — f^{uj). Here, 
F+g. (F~cr ) describes the tunneling rate of electrons with 
spin-cr into (out of) the QD from (into) the 77 lead with- 
out the occupation of the a state. Similarly, F+g. (F"^,.) 
describes the tunneling rate of electrons with spin-a into 
(out of) the QD, when the QD is already occupied by 
an electron with spin-tr, revealing the modification of the 
corresponding rates due to the Coulomb repulsion. The 
particle current /^ flowing from the lead 77 to the QD is 

Iilh = Z^i^i^aPdd + r,j^Pacr - F+^Pcrcr - F+„poo)- (3) 



This formula demonstrates that the current is totally de- 
termined by the diagonal elements of the density matrix 
of the central region. However, the nondiagonal element 
of the density matrix is coupled with the diagonal ele- 
ments in the rate equation (|2b|) . and therefore indirectly 
influences the tunneling current. 



III. QUANTUM SHOT NOISE FORMULA 

There is a well-established procedure, namely, the 
generation-recombination approach for multielectron 
channels, for the calculation of the noise power spec- 
trum based on the classical rate equations (classical shot 
noise) ,iSiiiiiS In this section we modify this approach in 
order to take into account the nondiagonal density matrix 
elements and derive the general expression for a quantum 
shot noise for the single QD. 

Before proceeding with investigation of current corre- 
lation, it is helpful to rewrite the quantum rate equations 
as matrix form: 



dp{t) 
dt 



Mp(t) 



(4) 



where p{t) = {pw-,P^\-,Pii-,Pdd, P\i, PuV is a vector 
whose components are the density matrix elements, and 
the 6x6 matrix M can be easily obtained from Eqs. ^. 
Correspondingly, we can write the average electrical cur- 
rents across the left [1^) and right [Ir) junctions at time 
t as: 



(/L(fl,)(0)--eE[^i(fi)PW 



(5) 



where fi, and f ^ are current operators and the 
summation goes over all vector [Fi(fl;)p(i)] elements 
(fc — 1,2,- •• ,6). The current operators contain the 
rates for tunneling across the left and right junctions 



espectivelyj^; and they can be read from Eq. J 


f, = ± 


/ r-^ r-^ \ 
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where the sign of the current is chosen to be positive 
when the direction of the current is from left to right, so 
that the + sign in the last equation is for rj — L and the 
— sign stands for rj = R. The stationary current can be 
obtained as: 



/ = e 



£[ 



^LP 



(0) 



^E[r«/ 



,(0) 
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0(-t)5][f,,f(-i)f,p(")]fe, (11) 



where 0{t) is the Heavisidc function and the two terms in 
Eq. ((TT|l stand for i > and for t < 0. The Fourier trans- 
form of propagator f{±t) is f{±Lo) = (^-1^1 - m\ 

where / is an unit matrix. We can further simplify this 
expression by using the spectral decomposition of the ma- 
trix M: 



M 



^ A„SE("")S-i = ^ APa 



(12) 



where p^"^ is the steady state solution of Eq. (0J and 



which can be obtained from 



Mp(°) = 0, 



(8) 



along with the normalization relation X]n=i Pn — 1- We 
would like to point out that in our quantum version of 
rate equations, it is easy to check ^^^ M„„j — (to — 
1,2,3,4), which implies that: 1) the Matrix M has a 
zero eigenvalue; 2) there is always a steady state solution 
p^"^; 3) the normalization relation J2n=iPni^) = 1 is 
independent on time. 

It is well known that the noise power spectra can be 
expressed as the Fourier transform of the current-current 
correlation function: 



Si^i^Alo) = 2 / dte-*[U.WV(0))-(/r,)(V)] 

J — OQ 

- 2(/^(t)V(0))^-2(/„)^(V)^. (9) 

A convenient way to evaluate the double-time correla- 
tion function is to define the propagator T{t) = exp[Mi], 
which governs the time evolution of the density matrix 
elements Pk{t)- The average value of the electrical cur- 
rents across the left (II) and the right {Iji) junctions at 
a time t is given by 



(/L(i.)(t)) = -e^[fi(;,)f(Op(°) 



(10) 



which allows us to switch the time evolution from the 
vector p{t) to the current operators. Thus, we identify 
i'L{R){t) = rL(i{)T(i) as the time-dependent current op- 
erators. With these time-dependent operators we can 
calculate correlation functions of two current operators 
taken at different moments in time. In particular, corre- 
lation function of the currents /,, and /,,' in the tunnel 
junctions 77 and 77', measured at the two times t and 
respectively, is given bjjiS 



(/,(t)V(o)) = 0(t)^[f,(t)f,,p 



(0)1 



where A are eigenvalues of the matrix M, S is a matrix 
whose columns are eigenvectors of M, E^""^ is a matrix 
that has 1 at nn place and all other elements are zeros, 
and Pa is a projector operator associated with the eigen- 
value A, so that r(±cj) is 



Pa 



f(±c.) = E T^ 



^iuj — X 



(13) 



Inserting expression for propagator T into Eq. I|ll|) 
current-current correlation in the w-space becomes 



(/,(i)V(0)).. = E 

A,fc 

+ E 



A,fc L 



A.fe 



iuj — A 



iV/\r^ (0) 



J k 



(14) 



Eventually, substituting Eq. H14|l into the noise definition 
Eq. © , and noting that summation over zero eigenvalue 
will be canceled out exactly by the term (/^)(/^'), we can 
obtain the final expression for a noise power spectrum: 



S'/„/„,(w) =^w"^r, 



Sch 



r„PAfvp("' 



f„,PAf„p'"' 



^-^ \ —iuj — A 

k,X^Q 



(15) 



where Sf."^^ is the frequency-independent Schottky noise 
originated from the self-correlation of a given tunneling 
event with itself, which the double-time correlation func- 
tion Eq. (|ll|l can not contain. Due to the fact that the 
current has no explicit dependence on the nondiagonal 
elements of the density matrix, it can be simply written 



as 



10 



s: 



Sch 



E|[f.p^°' 



(16) 



It can be shown that in the zero- frequency hniit, uj — 0, 
S{0) = Si^i^ = Si^i^ = -Si^i^ = -Sijii^. The Fano 
factor, which measures a deviation from the uncorrelated 
Poissonian noise, is defined as: 



F 



m 

2el ' 



(17) 



where 2e/ is the Poissonian noise. 



IV. NUMERICAL CALCULATIONS AND 
DISCUSSIONS 

In the foUowing we consider two magnetic configura- 
tions: the parallel (P), when the majority of electrons 
in both leads point in the same direction, chosen to be 
the electron spin-up state, a =t; and the antiparallel 
(AP), in which the magnetization of the right electrode 
is reversed. The ferromagnetism of the leads can be ac- 
counted for by means of polarization-dependent coupling 
constants. Thus, we set for P alignment 

Tlt - r^jt = (1 + p)r, Tli = Tr^ = (i - p)r, (is) 

while for AP configuration we choose 

Tlt = r^x = (1 + p)r, Tli = Tr^~ = (i - p)T. (19) 

Here, P denotes the tunneling coupling between the 
QD and the leads without any internal magnetization, 
whereas p {0 < p < 1) stands for the polarization 
strength of the leads. We work in the wide band limit, i.e. 
P is supposed to be a constant, and we use it as an energy 
unit in the rest of this paper. The zero of energy is chosen 
to be the Fermi level of the left and the right leads in the 
equilibrium condition (/i^ = fiji = 0). For clarity, the 
bias voltage, V, between the source and the drain is con- 
sidered to be applied symmetrically, /i^ — —fiR ~ eV/2. 
The shift of the discrete level due to the external bias is 
neglected. 

As a reference case for our analysis we use the analytic 
result for the case of paramagnetic electrodes, p — 0. 
This case is exactly solvable even for different couplings 
to the left and the right lead, Pl and P/j. The resulting 
Fano factor is spin- flip independent. In Coulomb block- 
ade regime, eV/2 < ed + U, it is given by 



F = 1 



4PLr/ 



(Pl + 2P/ 



(20) 



where the bias voltage is considered to be large (eV/2 ^ 
€d), SO that P^ == P+ == and P+ = Pl, P^ == Tr. The 
Fano factor depends only on the asymmetry in the cou- 
pling between the leads and the dot: it is equal to | for 
the completely symmetric case P^ = P^, and approaches 
to 1 when one of the coupling constants becomes much 
larger than the other one. 

In the opposite regime, when the energy e^ + C/ is far 
below the Fermi level /i^ {eV/2 ^ e^ -I- C/), we have 



P+ = 0, fz = f+ - and f+ = 
^ = V R, and the Fano factor is 


= r+ 


p - 1 TlPk 





(Pl + P/ 



(21) 



It is equal to ^ for completely symmetric couplings and 
to 1 for the asymmetric ones. 

Results of our numerical calculations for the current- 
voltage characteristic and the dependence on the Fano 
factor vs. bias voltage for P and AP configuration are 
presented in Figs. 1-3. In the following we set td = 1, 
Coulomb interaction U — A. By applying the bias voltage 
we are varying the fermi levels of the leads. Two steps in 
the current voltage characteristic occur: one is when the 
Fermi level of the source, hl , crosses the discrete levels 
Erf (for eV/2 > ed) and the other is when the Fermi level 
/ii crosses ed + U (for eV/2 > Cd + U). 

The effects of the polarization on the Fano factor with- 
out the spin-flip scattering {R ~ 0) are plotted in Fig. 1. 
An increase of the polarization will lead to an enhance- 
ment of the current noise in both configurations (P and 
AP) but for different reasons. Let us consider P and AP 
configurations separately. 
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FIG. 1: Current (a, b) and Fano factor (c, d) vs. the bias volt- 
age in the P configuration (a, c) and AP configuration (b, d) 
calculated for no spin-flip processes and different polarization 
p. Other parameters are: ed = 1, U — 4, T = 0.1. 

When the leads arc in a P configuration [Fig. l(a),(c)], 
an increase of the polarization increases the tunneling 
rates Pl| and P/j| for electrons with the spin-up and de- 
creases the tunneling rates P^j^ and Tri for spin-down 
electrons. This will increase the spin-up current and de- 
crease the spin-down current but it will not affect the 
total current through the system, which is equal to the 
sum of the spin-up and spin-down current. In the limit 



where the Coulomb interaction prevents a double occu- 
pancy of the dot, there will be competition between tun- 
neling processes for electrons with the spin-up and those 
with the spin-down. The characteristic time for these 
two processes, due to polarization, is unequal: there 
is fast tunneling of spin-up electrons and slow tunnel- 
ing of spin-down electrons through the system. The 
spin which tunnel with a lower rate modulate tunnel- 
ing of the other spin-direction (so-called dynamical spin- 
blockade) >iSiiLi& Eventually, for a large value of polariza- 
tion, it leads to an effective bunching of tunneling events 
and, consequently, to the supper-Poissonian shot noise. 

Increasing the bias voltage above the Coulomb block- 
ade regime, i.e, for eF^/2 > e^ + U, opens one more 
conducting channel and removes spin-blockade. In this 
regime, spin-up and spin-down electrons are tunneling 
through the different channels and there is no more com- 
petition between these two tunneling events. This leads 
to a reduction of the current fluctuation and the Fano 
factor becomes the same as in the paramagnetic case. 

The situation is completely different in the AP con- 
figuration [Fig. 1(b), (d)]. An increase of the polariza- 
tion increases tunneling rates Flj and Tri and decreases 
tunneling rates F^^ and Fflf . An electron with the spin- 
up, which has tunneled from the left electrode into the 
QD, remains there for a long time because the tunneling 
rate F/j| is reduced by the polarization. This decreases 
the spin- up current. An increase of the polarization also 
decreases the spin-down current because it reduces the 
probability for tunneling of the spin-down electrons into 
the QD. This will decrease a total current through the 
system. The enhancement of the noise in the AP config- 
uration is due to the asymmetry in the tunneling rates 
into and out of the QD {Tl] > ^ri but F^j < Tjii) for 
each spin separately. 

For large voltage, in the regime eV/2 > ed + U, both 
conducting channels become available which results in re- 
duction of the noise comparing with the Coulomb block- 
ade regime. In this case the Fano factor does not go to 
the paramagnetic value because the asymmetry in the 
tunneling rates are still present. 

The Fano factor dependence on spin-fiip scattering can 
be analyzed from Figs. 2 and 3. The spin-fiip scattering 
will open one more path in the tunneling of an electron 
out of the QD: an electron with a spin- up (down) can tun- 
nel into the QD, and due to the spin-fiip scattering, it can 
tunnel out of the QD as an electron with spin-down(up). 
An electron which was spending more time in QD due to 
polarization (spin-down in the P configuration and spin- 
up in the AP configuration) now has more probability to 
tunnel out from the dot. This causes the decrease of the 
current fiuctuation and Fano factor. However, this is not 
true for the P-configuration when double occupation of 
the QD is allowed (for eV/2 > e^ + U). In this regime 
the spin-flip does not have any effect on the Fano factor 
[Figs. 2(c) and 3(c)]. Spin-up and spin-down electrons 
are passing through separate channels without changing 
their spins. 
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FIG. 2: Current (a, b) and Fano factor (c, d) in the P con- 
figuration (a, c) and AP configuration (b, d) calculated for 
polarization p — 0.5 and different spin-ffip processes. Other 
parameters are the same as in Fig.0 
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FIG. 3: The same as Fig. |21except for p = 0.9. 



V. CONCLUSION 

In this paper we have analyzed the zero frequency cur- 
rent shot noise through a quantum dot connected to two 
ferromagnetic leads, taking account for the effects caused 
by the Coulomb interaction and weak intradot spin-flip 
scattering. For this purpose the rate equation approach 
for the description of the shot noise was modified to in- 
clude the coherent electron evolution inside of the dot. 



We calculated the Fano factor dependence on the lead 
configuration (parallel and antiparallel configurations), 
the degree of their polarization and the amplitude of the 
spin- flip process. Our results readily show the regimes 
when the information provided by shot noise measure- 
ment will significantly exceed the information which can 
be obtained from measuring the current. 
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